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NOMENCLATURE 


D  jet  pipe  diameter 

dp  particle  diameter 

r  radial  coordinate 

S  mean  Stokes  number  (t  /  t_  ) 
m  m  p 

Sg  mean  Stokes  number  at  jet  exit 

t  time  scale  of  mean  flow 

m 

tp  particle  time  constant 

U  streamwise  component  of  fluid  mean  velocity 
Uq  bulk  velocity  of  fluid  at  exit 

U  centerline  velocity  of  fluid  at  exit 
om 

Up  streamwise  component  of  particle  mean  velocity 

u  streamwise  component  of  fluid  rms  velocity 

u  streamwise  component  of  particle  rms  velocity 
P 

uv  cross-correlation  of  fluid  velocity 

u  v  cross-correlation  of  particle  velocity 
P  P 

v  radial  component  of  fluid  rms  velocity 
Vp  radial  component  of  particle  rms  velocity 

w  tangential  component  of  fluid  rms  velocity 

z  streamwise  coordinate 

v £  kinematic  viscosity  of  the  fluid 

Pp  particle  density 

p j  fluid  density 
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I.  INTRODUCTION 


Interior  ballistics,  and  especially  the  early  phase  of  the  ballistic  cycle,  is  marked  by 
the  prevalence  of  two-phase  flows.  As  yet  only  limited  progress  has  been  demonstrated 
in  the  measurement  of  spatially  resolved  velocity  distributions  of  the  phases  which  is  of 
interest  to  those  trying  to  model  and  understand  the  intricacies  of  the  interior  ballistic 
phenomena.  Specific  examples  include  primer  flow  characterization,  the  dynamics  of 
liquid  droplets  in  LP  (liquid  propellant)  gun  systems  and  the  effect  of  particulates  on 
the  erosion  of  gun  tube  walls. 

While  in-situ  measurements  under  actual  firing  conditions  with  the  instrumentation 
described  here  are  not  yet  possible,  building  block  type  studies  can  yield  considerable 
insight  into  the  nature  of  two-phase  flows.  Here,  we  report  on  such  an  approach, 
measurements  on  a  two-phase,  subsonic  jet  by  means  of  LDV  (Laser  Doppler 
Velocimetry).  The  mass  loading  was  1%  and  the  jet  exit  Reynolds  number  was  29  000. 
Velocities  of  both  phases  as  well  as  the  turbulence  quantities  were  measured  at  several 
locations  downstream  of  the  nozzle  exit  and  transverse  to  the  bulk  flow.  This  flow 
geometry  is  akin  to  the  situation  encountered  in  the  study  of  primer  flows  in  cartridge 
cases. 

The  report  is  divided  as  follows:  Section  II  presents  the  salient  features  of  the  flow 
and  the  parameters  used  to  characterize  it.  Next,  the  experimental  arrangement  and  the 
protocol  to  obtain  the  data  are  described.  This  is  followed  by  a  discussion  and 
interpretation  of  the  results  in  Section  IV.  The  last  section  draws  conclusions  from  this 
study. 


II.  ANALYTICAL  BACKGROUND 

Spurred  by  the  emergence  of  several  new  diagnostic  techniques,  including  LDV,  the 
study  of  turbulent,  dispersed  two-phase  flows  has  witnessed  a  resurgence  in  the  last  few 
years.  Several  groups  of  researchers  have  attempted  two-phase  flow  velocity 
measurements  in  particle  laden  flows.  The  most  notable  of  these  are  summarized  in 
Table  I  and  references  1-19.  All  of  the  studied  flows  are  characterized  by  relatively  low 
density  of  loading  and  moderate  Reynolds  number  and  narrow  range  of  particle  sizes. 
Only  in  a  few  cases  was  data  on  the  nature,  character  and  magnitude  of  the  turbulence 
presented.  Theoretical  approaches  are  discussed  in  Ref.  5,  20  and  21. 

■w 

The  theory  of  steady,  subsonic  jets  is  well  documented  and  allows  us  to  verify  the 
experimental  technique  and  the  reasonableness  of  the  results.  We  will  concentrate  on 
three  important  aspects  of  jet  behavior:  centerline  decay,  geometric  form  of  the  velocity 
profile  and,  with  the  second  phase  present,  the  value  of  the  Stokes  number,  a  measure 
of  the  responsiveness  of  the  particles  to  the  flow. 

Along  the  axis  of  the  jet,  the  velocity  decay  in  the  far  field  can  be  represented  by 


U/UQ  =  6.57  D/Z 

due  to  Spalding,  see  Ref.  11.  The  time  scale  of  the  flow  can  be  expressed  as 


(1) 


tm  =(dU/dZ  r1  =  (D/Uq)  (Z/D)2(l/6.57) 


(2) 
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where  again  U  is  the  centerline  velocity  at  the  point  of  observation  and  Uq  is  the  exit 
bulk  velocity  of  the  jet.  It  follows  that  the  time  scale  at  the  jet  exit  is  then  (D/Uq). 

In  addition,  a  particle  time  constant,  t  has  to  be  introduced.  Following  Snyder  and 
Lumley,  Ref.  21,  and  where  is  the  fluid  kinematic  viscosity, 


v = dP2(i+  2  y^t  >  / 36  "f  •  (3) 

With  these  parameters  a  mean  Stokes  number,  Sm  can  be  defined  as  follows: 

Sm  =  lm  /  ‘p  =  Se  l1/6'57)  <Z/D)2  ’  W 

where  Sg  is  the  mean  Stokes  number  at  the  jet  exit, 

Se  =  (D/U0)  /  tp  ,  (5) 

see  Hardalupas  et  al,  Ref.  11. 

The  presence  of  particles  in  the  flow  has  an  important  effect  on  the  development 
and  character  of  the  velocity  distribution.  The  degree  of  influence  of  the  second  phase 
can  he  quantified,  and  trends  predicted,  if  the  value  of  the  Stokes  number  is  known. 
Generally,  if  Sm  <  1,  the  particles  are  unresponsive  to  the  flow  and  a  mean  slip 

between  the  two  phases  occurs.  At  Sm  >  1,  the  particles  respond  to  the  flow.  In  general, 

the  momentum  flow  rate  of  the  gaseous  phase  is  affected  by  the  particulate  phase  only 
if  the  momentum  flow  rate  of  the  latter  exceeds  about  10%  that  of  the  carrier  phase.  In 
addition,  it  is  expected  that  particle-to-particle  interaction  will  take  place  for  volume 
loading  of  more  than  0.3%. 


III.  THE  EXPERIMENTAL  SETUP 

A  schematic  of  the  jet  flow  rig  with  the  instrumentation  used  for  the  measurements 
is  shown  in  Fig.  1.  The  jet  issues  vertically  downward  from  a  pipe  of  25.4  mm  in 
diameter  and  0.56  m  in  length,  giving  an  L/D  ratio  of  22,  to  a  plexiglass  enclosure 
which  eliminates  the  effects  of  room  air  currents  and  ensures  that  the  particles  are 
contained  and  safely  exhausted  out  of  the  room.  The  enclosure  is  relatively  large  with 
dimensions  of  16D  x  16D  x  16D  so  that  the  jet  flow  can  be  considered  to  be  unconfined. 
The  air  to  the  jet  was  supplied  from  a  fan  at  a  steady  rate  of  9.01  1/s  and  monitored  by 
a  flowmeter  within  1%  accuracy.  The  corresponding  Reynolds  number  at  the  jet  exit 
was  29  000  based  on  the  bulk  velocity  of  17.7  m/s.  The  particles  were  introduced  in  a 
plenum  upstream  of  the  jet  pipe  by  means  of  a  rotating  brush/piston  assembly,  TSI 
model  3410  particle  disperser.  The  particle  loading  was  controlled  by  the  piston  driven 
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variable  speed  motor  and  supplemented  by  the  air  flow  through  the  disperser  which  was 
provided  by  a  compressor  and  controlled  by  a  pressure  regulator.  For  the  single  phase 
flow  measurements  the  particle  disperser  was  used  to  seed  the  flow  with  titanium 
dioxide  particles  of  micron  size.  For  the  two-phase  flow  measurements  the  disperser 
delivered  glass  beads  of  80/im  nominal  mean  diameter  at  a  mass  loading  of  1%.  The 
glass  beads  had  a  relatively  narrow  size  range  of  60-95/im.  so  that  the  effects  associated 
with  polydispersed  particles  were  minimized. 

The  velocity  of  the  fluid  and  particle  flow  were  measured  by  a  two-component, 
dual-beam  laser-Doppler  velocimeter  mounted  on  a  computer  controlled  traversing 
mechanism  that  allowed  the  movement  of  the  measurement  volume  in  three  directions. 
The  LDV  consists  of  a  high  power  argon  laser  source  (2W)  and  optics  which  forms  two 
measuring  volumes.  The  backscattered  light  from  the  two  control  volumes  is  collected, 
color  separated  and  each  focussed  on  the  pinhole  of  a  photomultiplier.  The  Bragg  cells 
were  used  for  frequency  shifting  to  resolve  negative  velocities  in  both  components.  The 
principal  characteristics  of  the  laser  Doppler  velocimeter  are  summarized  in  Table  2. 
The  outputs  of  the  photomultipliers  were  band-pass  filtered,  amplified  and  input  to  two 
frequency  counters  which  were  interfaced  with  each  other  and  to  a  microcomputer  so  as 
to  allow  measurements  of  the  two  components  of  velocity  within  a  coincidence  time 
window  of  around  100/is.  The  microcomputer  stored  the  data  and  at  the  end  of  data 
acquisition  performed  statistical  calculations  to  yield  the  mean  and  rms  velocities  (U, 
V,  W  and  u,  v,  w  respectively),  and  the  cross-correlations  uv  and  uw. 

Consistency  of  the  data  was  repeatedly  verified  by  retaking  data  at  selected  points. 
The  major  sources  of  error  associated  with  the  velocity  measurements  are  due  to 
finite-size  statistics  and  velocity  gradient  broadening  effects,  Ref.  22.  The  overall  error 
in  the  mean  velocity  is  estimated  to  be  of  the  order  of  3%  and  around  10%  for  the  rms 
velocity  and  cross-correlation.  Of  course,  some  variation  in  the  radial  distribution  was 
noted,  especially  near  the  edge  of  the  jet,  due  to  the  changing  nature  of  the  flow.  The 
magnitude  of  these  errors  is  considered  to  be  too  small  to  affect  the  conclusions  drawn 
on  the  basis  of  results  presented  and  discussed  in  the  following  section.  Due  to  the 
sparsity  of  published  two-phase  flow  data  and  the  choice  of  flow  conditions  which  differ 
from  ours,  a  direct  comparison  with  the  work  of  others  is  not  possible. 


IV.  RESULTS  AND  DISCUSSION 


A.  Background 

In  this  section  sample  results  selected  to  depict  salient  features  of  the  flow  are 
presented  and  discussed  under  two  headings.  The  first  is  concerned  with  velocity 
characteristics  of  the  single-phase  flow  and  the  second  discusses  the  results  of  the  flow 
containing  the  particles  and  contrasts  these  with  the  single  phase  results.  The  velocity 
results  are  presented  in  the  form  of  centerline  distributions  and  radial  profiles  and  are 
all  normalized  by  the  mean  exit  centerline  velocity  U  of  around  20  m/s.  This  allows 
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not  only  a  clear  description  of  the  jet  centerline  velocity  decay  but  also  gives  a  better 
comparison  between  the  particle  and  single-phase  flow  velocities. 


The  centerline  velocity  distributions  of  the  streamwise  mean  and  rms  velocity,  the 
radial  rms  velocity  and  some  of  the  streamwise  velocity  probability  density  functions, 
also  called  pdfs,  obtained  at  various  downstream  locations  are  shown  in  Fig.  2.  The 
mean  velocity  is  constant  over  a  distance  of  4  jet  diameters  and  the  corresponding  rms 
velocity  also  remains  constant  over  this  distance  at  a  value  of  around  4%  of  the  mean. 
Further  downstream,  the  mean  velocity  decays  and  the  rms  velocity  gradually  increases 
to  around  19%  of  the  local  mean  at  z  =  9D.  The  radial  rms  velocity  also  increases  and 
is  always  less  than  the  streamwise  value,  as  expected,  by  about  30%.  The  velocity  pdf 
at  z  =  3D  is  near  Gaussian  with  a  relatively  small  standard  deviation  (rms),  but  at  z  = 
6D  the  pdf  becomes  skewed  due  to  the  intermittency  of  the  flow  as  the  shear  layers 
start  to  merge  together  at  around  z  =  4D.  At  farther  downstream  locations  the  two 
shear  layers  mix  more  readily  and  the  velocity  pdf  becomes  near  Gaussian  and  is  typical 
of  the  fully  developed  turbulent  jet  flow. 

A  more  detailed  description  of  the  flow  is  given  in  Fig.  3  which  shows  the  radial 
profiles  of  the  streamwise  mean  and  rms  velocity,  radial  rms  velocity  and  the 
corresponding  Reynolds  shear  stress  of  uv  at  three  downstream  locations  of  z  =  ID,  5D 
and  9D.  Fig.  3  is  consistent  with  Fig.  2  and  shows  that  the  mean  velocity  profile  close 
to  the  pipe  exit  (z  =  ID)  is  similar  to  fully  developed  turbulent  pipe  flow  with 
corresponding  rms  levels  of  around  4%  on  the  centerline  increasing  to  a  maximum  of 
around  14%  of  the  centerline  velocity  in  the  shear  layers.  As  expected,  the  radial  rms 
velocities  are  in  general  70%  of  the  streamwise  values.  The  profiles  show  good 
symmetry  about  the  centerline  at  all  three  downstream  locations.  The  velocities 
obtained  along  the  orthogonal  diameter  at  z  =  5D  also  indicate  an  excellent 
axisymmetry  of  the  flow;  the  transverse  rms  velocity  in  this  case  corresponds,  as 
discussed  above,  to  the  tangential  component  and  is  almost  identical  to  the  radial  rms 
velocity.  The  mean  velocity  profiles  become  flatter  at  downstream  locations  as  the  jet 
spreads  and  turbulence  tends  towards  homogeneity,  with  streamwise  rms  velocities  of 
around  20%  of  the  centerline  velocity  at  z  =  9D.  Parenthetically  it  should  be  noted  that 
most  of  the  published  data  on  single  phase  jets  deal  with  the  far  field,  so  direct 
comparison  with  the  data  reported  here  is  problematic.  However,  the  trends,  observed 
by  Wygnanski  and  Fiedler,  Ref.  23,  for  example,  on  centerline  velocity  decay  and  radial 
and  axial  rms  values  were  also  apparent  in  these  measurements. 


C.  Two-Phase  Flow 

The  velocity  measurements  of  the  particle  flow  are  presented  in  Fig.  4-6.  Here,  the 
corresponding  single-phase  flow  results  are  reproduced  to  allow  comparison  of  the  two 
flows.  It  is  assumed  for  the  relatively  dilute  concentrations  considered  here  that  the 
single-phase  flow  results  closely  approximate  that  of  the  carried  phase,  as  pointed  out  in 
Section  II.  The  results  of  Ref.  11,  for  example,  show  almost  identical  flow  fields  for  the 
single  and  carried  phase,  particularly  in  the  near  field  of  the  jet,  with  mass  loadings  of 
the  same  particles  up  to  20%. 
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The  Stokesian  time  constant  of  the  particles,  tp,  Equation  3,  is  about  50  ms  and  is 

much  higher  than  the  time  constant  of  the  mean  flow,  t  ,  Equation  2,  which  varies 

from  around  1.4  ms  at  the  exit  to  about  17  ms  at  z  =  9D;  the  mean  slip  between  the 
two  phases  is  expected  when  tm  <  tp,  i.e.  Sm  <  1.  Fig.  4  shows  the  centerline 

distributions  of  the  streamwise  mean  velocity  for  the  two  phases  and  indicates  that  the 
glass  beads  lag  the  fluid  flow  by  about  8%  near  the  exit,  caught  up  and  led  it  by  around 
7%  at  z  =  9D.  Fig.  4  also  presents  the  streamwise  and  radial  rms  velocities  together 
with  the  pdf's  of  particle  streamwise  velocity  at  various  downstream  locations. 
Although  the  radial  rms  velocity  of  the  beads  is  lower  than  that  of  the  fluid  at  the  same 
locations  along  the  centerline  indicating  a  relatively  small  cross-stream  spread  of  the 
beads,  the  streamwise  rms  values  are  surprisingly  higher,  up  to  3.5  times  at  the  exit, 
than  the  corresponding  values  of  the  fluid  in  the  region  down  to  about  z  =  7D.  The 
velocity  pdfs  are  consistent  with  this  observation  in  that  they  broaden  toward  lower 
velocities  compared  to  those  in  Fig.  2.  There  are  two  possible  explanations  for  this.  The 
first  one  may  be  due  to  the  polydispersed  nature  of  the  glass  beads,  albeit  they  are  in  a 
relatively  narrow  range  of  60-95/im.  As  a  consequence,  the  bigger  particles  with  the 
larger  time  constants  slip  more  and  move  more  slowly  than  the  small  ones.  The  second, 
and  probably  more  plausible  explanation,  is  that  of  Ref.  11,  and  referred  to  as 
"fan-spreading".  There  it  is  argued  that  the  trajectories  of  the  beads  at  the  exit  are 
virtually  straight  lines  rather  than  "rand,  m  walks"  since  in  their  passage  through  the 
pipe  leading  to  the  exit,  the  particles  bounce  off  the  wall  thereby  losing  most  of  their 
velocity  in  the  radial  direction.  The  angle  at  which  the  particles  spread,  outside  of  the 
pipe,  can  be  approximated  from  the  ratio  of  the  rms  velocity  in  the  radial  direction  to 
the  flow  mean  axial  velocity  which  typically  turns  out  to  be  around  ±3  degrees  in 
deviation  from  the  centerline.  In  a  plane  parallel  to  the  exit,  but  a  few  diameters 
downstream,  at  each  location  contribution  to  the  particle  motion  from  both  high 
velocity,  i.e.  particles  originally  near  the  centerline,  and  the  slower  particles  which 
originated  from  near  the  exit  wall,  can  be  expected,  leading  to  the  "broadening"  of  the 
velocity.  The  streamwise  velocity  at  downstream  locations  is  therefore  broadened 
towards  lower  velocities  by  contributions  from  beads  originating  near  the  wall, 
consequently  increasing  the  rms  values  as  observed  in  Fig.  4. 

Fig.  5  shows  the  comparison  of  the  streamwise  mean  velocity  profiles  for  the  two 
phases  at  three  downstream  locations.  The  bead  velocity  profiles  are  flatter  and, 
consistent  with  the  lower  radial  rms  values  of  Fig.  4,  indicate  a  lower  rate  of 
cross-stream  spread  of  the  beads.  The  corresponding  profiles  of  the  turbulence 
quantities  are  shown  in  Fig.  6.  As  in  Fig.  4,  the  streamwise  rms  values  at  the  first  two 
downstream  locations  are  higher  than  fluid  rms  velocities  near  the  centerline,  and 
everywhere  at  z  =  9D,  the  two  rms  values  become  comparable  and  support  the 
explanation  of  "fan-spreading".  The  cross-correlation  values  of  the  beads  (upvp)  are 

much  lower  than  the  fluid  values  as  was  also  noted  in  Ref.  11.  Thus  the  particles  at  this 
stage  are  only  subsidiary  contributors  to  the  total  turbulent  state  of  the  mixture. 


V.  CONCLUSIONS 


Measurements  in  an  unconfined,  axisymmetric,  80/r  particle  laden  jet  were  carried 
out  in  the  Interior  Ballistic  Flow  Diagnostic  Lab  using  LDA  techniques.  The  second 
phase  consisted  of  glass  beads  with  a  density  of  loading  of  1%.  For  comparison,  the  jet 
flow  parameters  were  also  measured  in  the  absence  of  the  second  phase. 
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The  most  important  findings  of  these  series  of  experiments  were  the  following: 

©  Single  phase  measurements  showed  excellent  axisymmetry  of  the  jet  and  were 
consistent  with  well  accepted  results  for  single  phase  jets. 

©  The  particles  having  a  Stokesian  time  constant  of  around  50  ms  were 

unresponsive  to  the  mean  and  turbulent  flow  and  lagged  the  fluid  by  about  8% 
near  the  jet  exit  but  later  led  it  by  about  7%  at  9  jet  diameters  downstream. 

©  The  mean  velocity  profiles  of  the  particles  were  flatter  than  those  of  the  fluid 
and  together  with  the  lower  values  of  the  radial  rms  velocity,  indicate  a 
relatively  small  cross-stream  spread  of  the  beads. 

©  Surprisingly,  the  streamwise  velocity  fluctuations  of  the  particles  near  the 
centerline  and  exit  were  much  higher,  up  to  3.5  times  at  the  exit,  than  the 
corresponding  fluid  values  in  a  region  of  up  to  7  jet  diameters.  This 
observation  is  consistent  with  the  findings  of  other  investigators  and  is  mainly 
attributed  to  the  fan-like  spreading  of  the  beads  in  a  straight-line  trajectory 
from  the  jet  exit. 

More  experiments  encompassing  variations  in  jet  exit  diameters,  density  of  loading 
and  the  interaction  of  particle  laden  jets  are  under  way  and  will  be  reported  in  the  near 
future. 
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Table  I.  REPRESENTATIVE  TWO-PHASE  FLOW  MEASUREMENTS 


Reference 

Flow  Configuration 

Particle  and 
Particle  Size 

Reynolds 

Number 

Particle  Mass 
Fraction 

Pepper,  Abuaf, 
Hetsroni  (1974) 

round  jet 

oil,  50 /i 

10  600 

49  300 

1  x  10~3 

Levy  and  Lockwood 
(1981) 

round  jet 
x/D  =  20 

sand,  215-1060/1 

25  000 

1.14-3.5 

Tsuji,  Morikawa 
(1982) 

pipe  flow 

30  mm  i.d. 

200/i,  3400 /i 

5  000 

not  given 

Wells  and  Stock 
(1983) 

square  wind 
tunnel 

glass  beads 

5/i,  57/i 

0.887 

(particle) 

n.a. 

Shuen,  Solomon, 
Zhang,  Faeth 
(1984) 

round  jet 
x/D  <  50 

sand,  79-207/1 

19  000 

0.2-0.66 

Modaress,  Tan, 
Elghobashi  (1984) 

round  jet 

glass  beads 

13  300 

0.32,  0.85 

Vames,  Hanratty 
(1988) 

round  jet 

water  droplets 
50-150/1 

96  000 

not  given 

Tsuji  et  al 
(1988) 

round  jet 

polystyrene 

170-1400M 

32  000 

not  given 

present  work 
(1989) 

round  jet 
x/D  =  20 

glass  beads, 

80/i 

29  000 

1.0 
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Table  II.  CHARACTERISTICS  OF  THE  LASER-DOPPLER  ANEMOMETER 


Laser  wavelength  514.5nm 

Fringe  spacing  4.5  /an 

Length  of  control  volume  at 
2 

1  /  e  intensity  0.97  mm 

Diameter  of  control  volume  at 
1  /  e  intensity  56.0  /an 

Number  of  fringes  at 

1  /  e^  intensity  12 

Frequency  shift  1-5  MHz 


Table  III.  FLOW  CONDITIONS  AT  JET  EXIT,  Z  =  1  D 


Jet  diameter 

2.54  cm 

Carrier  phase  exit  velocity  (peak) 

20.0  m/s 

Carrier  phase  exit  velocity  (bulk) 

17.7  m/s 

Turbulence  intensity 

1.4% 

488  nm 
4.3  /an 

0.93  mm 

53  /an 

12 

5  MHz 
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AIR  FLOW 
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FIGURE  1.  LAYOUT  OF  THE  EXPERIMENTAL  APPARATUS 


Figure  2.  Centerline  distribution  of  mean  and  rms  velocities  in  single-phase  flow. 
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Figure  3.  Radial  distributions  of  streamwise  mean  and  rms  velocities,  radial  and 
tangential  rms  velocity  and  cross-correlation  of  single-phase  flow. 
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Figure  5.  Radial  distributions  i,f  streamwise  mean  velocity  of  particle  flow. 
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Figure  6.  Radial  distributions  of  streamwise  rms  velocity  and  cross-correlation  of 
particle  flow. 
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